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ABSTRACT 

We estimate the solar system motion relative to the cosmic microwave background 
using type la supernovae (SNe) measurements. We take into account the correlations 
in the error bars of the SNe measurements arising from correlated peculiar velocities. 
Without accounting for correlations in the peculiar velocities, the SNe data we use 
appear to detect the peculiar velocity of the solar system at about the 3.5 a level. 
However, when the correlations are correctly accounted for, the SNe data only detects 
the solar system peculiar velocity at about the 2.5 a level. We forecast that the solar 
system peculiar velocity will be detected at the 9 a level by GAIA and the 11 cr level 
by the LSST. For these surveys we find the correlations are much less important as 
most of the signal comes from higher redshifts where the number density of SNe is 
insufficient for the correlations to be important. 



INTRODUCTION 



The cosmic microwave background (CMB) has a 3.4 
mK dipole anisotropy ((Hinshaw et al. ,2007.) which can 
naturally be explained as being due to the motion of 
the solar system with respect to the CM B rest frame 
llLvnd cn-Bcll ot al. 1989; Strauss ct al. 1992; Erdogdu et all 
I2OO6I : iLoeb fc Naravanl |2007| ). An interesting consistency 
check of this is to evaluate the solar sy stem motion from 
peculi ar velocity surveys (see for example iDale fc Giovanellil 

SNe luminosity measurements provide an accurate 
probe of peculiar velocities. Using observed correlations be- 
tween SNe light curves, we can estimate the SNe absolute 
magnitudes and thus obtain accurate distance estimates to 
the SNe. Combined with spectroscopic measurements of the 
host galaxies' redshifts, this can be used to estimate the pe- 
culiar velocity of each SNe's host galaxy. The motion of the 
solar system will then show up as a dipole anisotropy in the 
SNe derived peculiar velocities. It is interesting to compare 
the estimates of the solar system motion from the SNe with 
those derived from the CMB. If they turn out to be inconsis- 
tent then it may be an indication that there is a significantly 
large intrinsic temperature dipole on the CMB surface of 
last scattering (Turner 1991; Langlois & Piran 1996), which 
could be caused by a double inflation model l|Langloislll996l ) 
for example. 

A number of studies have made this comp arison 
ijRiess et al.lll995l : iBonvin et alllioo^ : Ijha et al.ll2007h . and 
a simplifying assumption used in these studies was that 
the peculiar velocities of the individual SNe were uncorre- 
lated with each other. However, as the peculiar velocities are 
caused by variations in the density field, neighbouring SNe 
will have correlated peculiar velocities (|Wang et al.l Il99^ : 



Sueiura et al. 1999'; iHui fc Green3l2006l : iBonvin et aDl2006l : 



Gordon et al.,i2007j ). These correlations will increase the er- 



ror bars on our peculiar velocity estimate as each new SNe 
measurement does not represent a completely independent 
realization of the velocity field. 

In this article we include the correlations of the peculiar 
velocities when estimating the motion of the solar system 
with respect to the cosmic rest frame. In Sec. [2] we give a 
simple example of the underestimation of the uncertainty 
that occurs when correlations between observations are not 
taken into account. In Sec.[3]we outline the formalism we use 
for the SNe correlations, and in Sec.|4]we apply the method 
to SALT calibrated SNe data. In Sec. [5] we look at the im- 
plications for future surveys. A summary and discussion of 
the results is given in Sec. EJ 



2 SIMPLE EXAMPLE OF CORRELATED 
ERRORS 

In order to illustrate the effect of correlated errors we con- 
sider a simple example (also discu s sed in Eq. 5 and its below 
parag raph of ICoorav fc Caldwelll (|2006l )) where we analyse 
data points (xi) drawn from a multivariate Gaussian like- 
lihood 



C (X \C\ "'"''^ exp(— (x — /x)'^C """(x — /i)/2). 



(1) 



The vector x is made up of the data points (xi) and each 
element of the vector /x is equal to a constant, /i. The co- 
variance matrix (C) has diagonal terms which are and the 
off-diagonal terms which are pa^. That is each data point 
has correlation p with the other data points. Suppose that 
one attempts to estimate the value of p and from the 
data and ignores correlations (i.e. assume p to be zero). The 
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maximum likelihood estimators of the mean and variance 
are then respectively given by 



(2) 



As the likelihood is a multivariate Gaussian distribution 

{xi) = /i and {{xi ~ fJ.){xj — fi)) = dj (3) 

where angular brackets denote the expectation value. Eval- 
uating Eq. ([2]) using Eq. (|3]) gives 

2 



{m) = fi, {s"^) ^ [1 - p)a^ 



(4) 



in the large N limit, where the approximation approaches 
equality for large A*'. As can be seen from the above equation, 
m is an unbiased estimator of the true mean, even when 
there are correlations in the data. However, as can also be 
seen, is biased by a factor of (1 — p). 

If there where no correlations, the error on jj, can be 
estimated by 

= ^ • (5) 



N 

If there are correlations present then it follows from Eq. Q 
that this estimator will give 



N 



(6) 



Now we find the correct value of Sjj. to see the effects of 
ignoring correlations. For large N, the expectation of the 
one sigma error on /x can be approximated 



5m ^ V(F-i)m,m (7) 

where {f~^)fi,^ is the {fJ,,ii) component of the inverse of the 
Fisher matrix. For a multivariate Gaussian likelihood the 
Fisher matrix is given by (see for example iXegmark et al.l 



dpadp0 / 



(8) 



(9) 



where a and P run over the different model parameters 
which are being estimated (fi and a, with p assumed known, 
in the current example). Eq. (O gives the true error on the 
value of p., and in our example this is 



1 + {N-I)p 



N 



(10) 



Comparing Eqs. ((6} and (|10|) we see that if the data are 
correlated (p > 0) but correlations are neglected then one 
would underestimate the uncertainty on p. One would over- 
estimate the error if the data is anti-correlated, although we 
note that this case is restricted as for the covariance matrix 

to be positive definite, p > —{N — 

As was shown in a earlier studie s (|Hui fc Greene! [200^ : 
iNeiU et al.|[2007l : [Gordon et aLll200'i1 ). an analogous under- 
estimation of the error happens if the correlations in low 
redshift SNe are not accounted for when using them in a 
sample to estimate the dark energy equation of state, w. 
In this article we show that there is also an underestima- 
tion in error on the motion of our solar system when the 
correlations in the SNe are not accounted for. 



3 METHOD 

The luminosity distance, (11, to a SN at redshift z, is defined 
such that 



C 



where is the observed flux and £ is the SN's intrinsic 
luminosity. Astronomers use magnitudes, which are related 
to the luminosity distance (in megaparsec) by 

jobs 



m — M = 5 logj^Q -I- 25, 



(11) 



where m and M are the apparent and absolute magnitudes 
respectively. In the context of SNe, M is a "nuisance param- 
eter" which is completely degenerate with log(i/o) and is 
marginalised over. For a Friedmann-Robertson- Walker Uni- 
verse the predicted luminosity distance is given by 

rf4.)Ml + ^)f ^ (12) 

(taking c = 1), where H is the Hubble parameter. In the 
limit of low redshift this reduces to dz, « z/ Ho. 

Given very stringent limits on the curvature of the uni- 
verse, we can safely work within the assumption of a flat- 
ness as the allowed curvature would not play any role at 
the scales of interest. In this case, the effect of a pecu- 
liar velocity (PV) leads to a perturbat i on in the luminosity 
distance (5dL) given by ('Sasaki"l987': 'Sugiura ct al .lll999|: 
Pvne fc Birk inshaw 2004: Bonvi n et al.,2006, : .Hui fc Greene! 
2006! ) 



Ml 



(l+ill, 

H{z) dL ' 



(13) 



where r is the position of the SN, and vo and v are the pecu- 
liar velocites of the observer and SN repectively. In the limit 
of low redshift, 5dL ~ r • [vo — v]///o- This demonstrates 
how a SNe survey that measures m and z can estimate the 
projected PV field. We now relate this to the cosmology. 

The projected velocity correlation function, 5(r,r') = 
((v(r) ■ f)(v(r') • f')), must be rotationally invariant, and 
therefore it can be decornposcd into a parall el and per- 
pendicular com ponents (jGorski, ,1988, : .Groth et al.l 1 19891 : 
lDodelson![2003 ) : 

^(r, r ) = sin 6 sin 0'C±( Ar, z, z')+cos 6 cos ei'^y (Ar, z, z) (14) 

where Ar = r — r', Ar = |Ar|, cos 9 = f • Af , and cos 9' = 
r' ■ Af . In linear theory, these ar e given by (|Gorski|[l988l : 
iGroth et al1ll989l : lDodelsonll2003! ): 



^Il,x=i?'(^) D'iz') 



poo 



^P(fc)i^ll,x(fe 



;r) 



(15) 



where for an arbitrary variable x, Kn{x) = joix) — 
K±{x) = ji{x)/x. D{z) is the growth function, and deriva- 
tives are with respect to conformal time. P{k) is the matter 
power spect rum which can be evaluated either numerically 
(e.g. CA MBlLewis et all bOO(t)) or using analytical approx- 
imations (|Eisenstein fc Hu!ll998l '). 

The above estimate of ^(r, r') is based on linear the- 
ory. On scales smaller than about 10/i~^Mpc, nonlinear con- 
tributions dominate. These are usually modeled as an un- 
correlated term which is independent of redshift, often set 
to (j„ ~ 300 km/s. Comparison with N-body simulations 
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ISilberman et al.ll200ll ) indicate that this is an effective way 
of accounting for the non-linearities. Other random errors 
that are usually considered are those from the lightcurve fit- 
ting (merr), aud intrinsic magnitude scatter {a,n) . It is just 
these three errors that are usually included in the analysis 
of SNe. 

The residual deviations of luminosity distance from the 
homogeneous expansion can be packed into a data vector 



'SdL\ ^ dt'ii)^ dL{z{i)) 



whose covariance matrix (from the correlated PVs) is given 

by 



(1 + ^) 
HdL 



(1 



HdL 



e(ri,rj), (17) 



while the standard uncorrelated random errors are given by 



( g I (o-.m+mcrr(l) )+( 1 



HdL 



(t;.(i8) 



Some example plots of Cv were given in iGordon et al.l 
l|2007l) . The hkeUhood is then 



where 



-/^exp(-iA^r-iA) 



T{i,j) = Cv{i,j) +cr{i) Si^ 



and 



(5dL\ 
V dL J, 



(1 + ^)^ 

H(z) di 



Ti ■ VQ. 



(19) 



(20) 



(21) 



We now proceed to find constraints on the observer ve- 
locity vq. We assume a standard ACDM cosmology and 
impose Big Bang Nucleosynthesis (BBN) prior Q},h^ ~ 
A/'(0.0214, 0.002) (Kirkman ct al. 2003), an d a Hubble Space 
Teles cope (HST) prior h ~ A/'(0.72, 0.08) jFreedman et all 
I2OOII ). These two priors remove models that are wildly at 
odds with standard cosmological probes, but do not unduly 
bias results towards standard cosmology. The likelihood has 
almost negligible dependence on , and to keep it in a range 
consistent with CMB and large scale structure estimates we 
give it a uniform prior n € [—0.9, 1.1]. 

We parameterize the solar system peculiar velocity as a 
magnitude (vo) and direction in galactic coordinates {l,b). 
The prior on {I, b) was assumed to be uniform on the sphere, 
i.e flat on I and cos (6 ). The prior on vq was set to be uniform. 

We use a SALT (|Guv et al.ll2005^ cahbrated low redshift 
SNe data selQ with heliocentric redshifts in the range cz G 
[2278,37163]km/s, {z G [0.0076,0.124]). A histogram of the 
61 redshifts is shown in Fig. [T]and the sky positions of the 
data are shown in Fig. (2] We also used the higher redshift 
71 SNe from the SNLS data set0 

The SALT calibration involves the additional param- 
eters (a, 13), which account for the shape/luminosity and 
colour/luminosity relations of SNe. These and the other 
parameters o"8, o"m, M) are all given broad uni- 

form priors. We use the standard Markov Chain Monte 
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Figure 1. Heliocentric distribution of redshifts for the low z sam- 
ple. 




Figure 2. Positions of the SNe for the low z sample in galactic 
coordinates. The size of the disk is inversely proportional to the 
redshift. The color of the disk is related to the relative luminosity 
distance error {ScLl / cLl)- The stars indicate the direction of the 
CMB dipole. 



http: / / snls.in2p3.fr/conf/ papers/cosmol/| ~ 



Carlo (MCMC) method to generate samples from the pos - 
terior distribution of the parameters (|Lewis fc Bridle|[200^ ). 
Convergence was checked using multiple chains with dif- 
ferent starting position s, and also the R — 1 statistic 
l|Gelman fc RubinI Il992l l. We also checked that the esti- 
mated posterior distributions reduced to the prior distribu- 
tions when no data was used. The analysis was checked with 
two completely independent codes and MCMC chains. 

Additionally, we looked at the combination of the 
SNe observations with the WMAP data of the CMB 
(Spergel et al.. ,2006 ) (with the usual CMB priors in this 
case) . We stress that we do not use the WMAP dipole infor- 
mation, but rather just the I > 1 information, which when 
combined with the SNe has the effect of constraining matter 
density and the amplitude of matter fluctuations. 



4 Gordon, Land, and Slosar 



1000 



>° 500 




Figure 3. Marginalized one and two sigma contours for the mag- 
nitude of the solar system peculiar velocity {vq ) vs the dispersion 
of the matter density field smoothed on scales of 8/i~^Mpc (erg). 
The thick (black) contours are for SNe data with BBN and HST 
priors. The thin (blue) contours are for SNe with a WMAP prior. 



4 RESULTS 

In Fig. [3] we plot the marginalized probability contours for 
vo and (Tg, where we see that higher values of erg have 
broader contours on vq- This is because a larger ag implies 
more correlations between the SNe peculiar velocities and 
so less of a reduction in the errors due to averaging effects. 
The contours for when WMAP is included are also plotted. 

In Fig.|4]a marginalized probability distribution is plot- 
ted for Vo, and we see that the uncertainty on vq increases 
dramatically when the correlations are accounted for. It can 
also be seen that adding WMAP data has a negligible effect. 
As seen from Fig.[3l this arises because the WMAP data con- 
strains (Tg « 0.8, which happens to lie on an approximately 
average value for the vq uncertainty. In Fig.[5]the one sigma 
confidence intervals are plotted for the direction of the solar 
system motion. As can be seen, not accounting for the cor- 
relations underestimates the uncertainty by about a factor 
of 2. 

Also, Figs. |4] and [5] shows that the SNe data are con- 
sistent with the CMB dipole estimat e of (vp, I, b) = (369 ± 
3km/s, 263.86 ± 0.04°, 48.24 ± 0.10°) (|Hinshaw et al.ll2007l ). 
In Table [1] the mean and uncertainties of the solar system 
peculiar velocity are given. We find that when the corre- 
lations are not included, the estimate of vq is about 3.5 
standard deviations from zero. While if the correlations are 
accounted for then vq is only about 2.5 standard deviations 
from zero. One can convert the number of standard devi- 
ations of th e detection into upper b ounds on the Bayesian 
odds ratio (jCordon fc Trottal boOTl ). Without correlations 
the odds, from SNe data, of vo being non-zero appear to 
be at best 119:1. While if the correlations are accounted for 
then the odds are at best only 7:1. In Table [1] we also give 



200 400 600 800 

Figure 4. Marginalized probability distributions for the magni- 
tude of the solar system peculiar velocity, vq ■ The dotted line is 
for when the correlations in the SNe peculiar velocities are not 
accounted for. The solid line is with correlations and the dashed 
line is with correlations and WMAP temperature data included 
to reduce the uncertainty in the cosmological parameters. 




Figure 5. One sigma contours for the direction of the solar 
system velocity. The cases plotted are when correlated (solid) 
and uncorrelated (dashed) SNe peculiar velocities are used. The 
star shows the direction as determined by the CMB. 



an estimate for the local group motion which was obtained 
by sub tracting the solar system velocity relative to the local 
group (|Yahil et al.lll977h . 



5 FORECASTS 

In this section we forecast constraints on the motion 
of the solar system from GAIA (the "super-Hipparcos" 
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Solar System uncorrelated 
Solar System correlated 
Local Group uncorrelated 
Local Group correlated 



I 

238 ± 26° 
234 ± 44° 
260 ± 14° 
257 ± 24° 



45 ± 14° 
39 ± 21° 
32 ± 11° 
29 ± 16° 



VQ (km/s) 

475 ± 134 
468 ± 186 
697 ± 137 
690 ± 201 



Table 1. The mean and standard deviation for the estimate of 
the solar system and local group velocity from current SNe data. 
The results for both the correlated and uncorrelated peculiar ve- 
locities are shown. 



VQ (km/s) 



GAIA Uncorrelated PVs 


8° 


5° 


36 


LSST Uncorrelated PVs 


7° 


5° 


32 


GAIA Correlated PVS 


10° 


6° 


42 


LSST Correlated PVS 


8° 


5° 


34 



Table 2. The forecasted marginalised standard deviation for the 
estimate of the solar system peculiar velocity from the future SNe 
surveys GAIA and LSST, where SNe peculiar velocities (PVs) are 
treated as uncorrelated and correlated. 



number 



1000- 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 



Figure 6. A plot of the lower bound on the number of SNe needed 
at or with redshift less than z for correlations in the peculiar 
velocities to be important (solid line). Also plotted are the number 
of SNe, with a redshift less than or equal to z, for the currently 
available low redshift SNe (dotted line) and for volume weighted 
GAIA (short dashed) and LSST (long dashed) surveys. 



satellite) a nd LSST. For GA I A, ba sed on the simu- 
lations by iBelokurov fc Evand (|2003l ). we generate a 
sample of 6,317 SNe distributed over the full sky with 
z < 0.14. For LSST we generate 30,000 SNe distri buted 
over the full sky with z < 0.3 (jWang et all l2005l l. We 
weighted the distribution of SNe by cos{b)z'^ to ac- 
count for the volume in spherical coordinates, that is 
we keep the density constant with z. As our fiducial 
model we took {vo, I, b, Qm, i^b, h, 

{369, 264°, 48°, 0.3, 0.041, 0.72, 0.96, -1, 0.85, 300, 0.1, 0.1}. 
We do not include the SALT calibration parameters (a, jS) 
in the forecast, but we have checked using forecasts for the 
data sets used in Sec. 3 that this does not have a significant 
effect. In order to use the Fisher matrix (see Eq. (|^), we 
consider the function 

d = df^ 10*'°/^ (22) 

where d°^'^ is given by Eq. (|lip . The expectation value vector 
has each element given by 

{d) = 10"«/^ {df + f • vo ^^^) . (23) 

with given by Ea. (ll2|l . The covariance matrix 

(C) is as before, Ea. (|20|l . but with the extra factor 
{iQ^'o/^Y Note the reason for the slightly dif- 

ferent function of the data compared to Eq. (|2ip is so that 
the data vector does not depend on any of the parameters. 

In Table [2] we present our main forecast results. As can 
be seen there is a dramatic improvement in the constraints 
compared to current data. Also, unlike current data, taking 
into account the correlations does not have a large effect. 
This is because most of the constraining power for GAIA 
and LSST comes from higher redshifts, where the peculiar 
velocity errors are negligible compared to the other types 
of error. This can be understood as follows. From Eq. [T7] 
we see that the error induced by peculiar coherent veloc- 
ity fiows drops as 1/z, and thus for high enough redshift 
(z > 0.015 for a typical experiment) they become unim- 
portant compared to the redshift independent measurement 
errors in Eq. (|18|) . In this limit the weakening of the dipole 



signal, that drops as 1/ z in Ea. (|2ip . is exactly compensated 
by the number of supernova in a redshift slice, which in- 
creases as z'^, for a volume- weighted survey. The signal to 
noise is therefore low and increasing at low redshift, tailing 
off to a constant value at redshifts at which peculiar veloci- 
ties become unimportant. 

As a rough measure of whether the correlated error will 
be important, at a redshift z, we look at the ratio N^a'^/Ci, 
where Ci, is evaluated using Eq. (|17|l with two SNe both at 
redshift z and 90° apart. This is effectively the ratio of the 
measurement error on the SN luminosity to its covariance 
with a typical SN in the dataset. Eq. (|10|l tells us that this is 
approximately equal to the number of SNe for which the co- 
variance will become important to our error estimates, and 
we plot this in Fig. [6] As can be seen this simple estimate is 
in agreement with the more complete analysis: current data 
are effected by correlated errors much more so than GAIA 
and LSST, which will be practically uneffected by correla- 
tions. In other words, under the assumption of a volume 
weighted redshift distribution for the SNe, GAIA and LSST 
can rely on the much higher redshift data to constrain our 
peculiar motion, as these are considerably less affected by 
PV correlations. 



6 DISCUSSION 

To summarize, we have used SALT calibrated SNe data to 
estimate the motion of he solar system. As seen from Ta- 
ble [1] the error bars are under-estimated by about 50% if 
the correlations in the peculiar velocity are not accounted 
for. 

We now compare our fin dings to previous published re- 
sults. In lBonvin et al.l (120061 ) they used 44 SALT cahbrated 
SNe. They only allowed {l,b,vo} to vary and all the other 
parameters where fixed to standard values. They did not ac- 
count for correlations in the peculiar velocities. They found 
Vo — 40 5 ± 192km / s whic h is compatible with our result. 

In IJha et all (|2007l ) they used 69 SNe with z € 
{0.005,0.025}. They also only allowed {l,b,vo} to vary. 
Additionally, they used MLCS2k2 to calibrate the data. 
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rather than the SALT method. They also did not account 
for the correlations in the peculiar velocities. They evalu- 
ated the motion of the local group and found {l,h,vo} — 
{258 ± 18°, 51 ± 12°, 541 ± 75km/s}. Our local group veloc- 
ity resul ts, which a re list ed in Table [1] are compatible with 
those of iJha et all (|2007l ') but, even when we don't take into 
account the correlations in the peculiar velocities, our error 
on t he magnitude of vo is about 80% larger than that of 
|jha et al.. (,2007) study. This is due to several factors. They 
had a lower redshift limit than us: 0.005 vs. 0.0076. One of 
the reasons we did not go to such a low redshift is that the 
peculiar velocities (including the motion of our solar system) 
become of order the Hubble expansion. This means that the 
motion of our solar syst em has a h i gh sig nal (hence the low 
error bars obtained by |jha et al.l i|2007l) ) but one can no 
longer use Eq. (|13|l to evaluate the effects of peculiar veloc- 
ity on the luminosity distance. It would be possible to use 
a higher order v ersion of Eq. (|13p but this was not done by 
IJha et al.l (|2007l ') and so their results will have an additional 
unreported systematic error due to the induced luminosity 
change being calculated incorrectly. Also, the extra very low 
redshift SNe that were used are dominated by SNe that are 
too close together to model the correlations in the peculiar 
velocity using linear theory, Eq. I\15\i . Overall, our results 
are the only ones that take into account the correlations in 
peculiar velocities and account for the uncertainties in the 
cosmological and calibration parameters. 

We also made forecasts for the GAIA and LSST sur- 
veys, assuming volume weighting for the redshift distribu- 
tion. We found that the error bars will be about 4 times 
smaller than those of current data, but still not competitive 
with those from the CMB by a factor of ~ 10 (assuming 
that the CMB dipole is due to our local motion). Also, for 
GAIA and especially LSST we found that correlations had 
little effect as most of the signal came from higher redshifts 
where the correlations are almost negligible for the sample 
sizes considered. 

In future work, we plan to test the techniques we have 
used in this paper against simulated SNe surveys generated 
from N-body simulations. These will test the assumptions 
that go into our d ata modeling - mos t importantly the effect 
of non- Unearities ( Haugboelle et al.|[2006i ') . 

In IWatkins fc FeldmanI (|2007D it was shown that the 
large sca le properties (bulk flow and shear (jKaisej Il99ll : 
Ijaffe fc Kaiser 199^)) of the peculiar velocity field derived 
from a low redshift sample of 73 SNe was consistent with 
the bulk flow and shear of the velocity field derived from 
the SFI, ENEAR, and SBF surveys. It would be interesting 
to combine all these surveys t ogether (with SFI replaced by 
SFI++ ((Springob et al.ll2007l )) to estimate the solar system 
velocity with respect to the CMB and put robust limits on 
the intrinsic CMB dipole. 
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